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ABSTRACT 

The measurement of redshifts for Gamma-Ray Bursts (GRBs) is an impor- 
tant issue for the study of the high redshift universe and cosmology. We are 
constructing a program to estimate the redshifts for GRBs from the original 
Swift light curves and spectra, aiming to get redshifts for the Swift bursts with- 
out spectroscopic or photometric redshifts. We derive the luminosity indicators 
from the light curves and spectra of each burst, including the lag time between 
low and high photon energy light curves, the variability of the light curve, the 
peak energy of the spectrum, the number of peaks in the light curve, and the 
minimum rise time of the peaks. These luminosity indicators can each be related 
directly to the luminosity, and we combine their independent luminosities into 
one weighted average. Then with our combined luminosity value, the observed 
burst peak brightness, and the concordance redshift-distance relation, we can de- 
rive the redshift for each burst. In this paper, we test the accuracy of our method 
on 107 bursts with known spectroscopic redshift. The reduced of our best red- 
shifts (zbest) compared with known spectroscopic redshifts (zspec) is 0.86, and the 
average value of logio^Zhest/ Zgpec) is 0.01, with this indicating that our error bars 
are good and our estimates are not biased. The RMS scatter of logiQ{z},est/ ^spec) 
is 0.26, with a comparison of 0.30 for RMS of logiQ{zspec)- We made a selection on 
bursts with relatively accurate redshift estimation. The RMS of log {zbest / z spec) 
decreases to 0.19, and the RMS scatter of logio{zspec) for this subsample is 0.28. 
For Swift bursts measured over a relatively narrow energy band, the uncertainty 
in determining the peak energy is one of the main restrictions on our accuracy. 
Although the accuracy of our Zbest values are not as good as that of spectroscopic 
redshifts, it is very useful for demographic studies, as our sample is nearly com- 
plete and the redshifts do not have the severe selection effects associated with 
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optical spectroscopy. 

Subject headings: Gamma-Ray: Bursts - Cosmology: observations 
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1. Introduction 

The redshifts have long been an important issue for long-duration Gamma-Ray Bursts 
(GRBs). Since the first X-ray, optical, and radio counterparts were discovered in 1997 
(Costa et al. 1997; van Paradijs et al. 1997; Frail et al. 1997), GRBs are confirmed to be in 
galaxies at cosmological distances. However, until now, only a small fraction of the bursts 
have their redshifts measured. Even Swift, over its first four years of operation, has only 
roughly 30% of its bursts with spectroscopic/photometric redshifts. These redshifts may 
have complex selection biases as a function of redshift relating to the difficulty of getting 
redshifts for faint and distant bursts as well as the distribution of bands used for spectra 
and photometry. One illustration of the importance of the selection effects is that the 
average of redshift value for Swift GRBs have been dropped from < z > = 2.8 (Jakobsson 
et al. 2006c) and to < 2; > = 2.1 (Jakobsson 2008) while the average redshift value of 
earher sateUites is much lower yet (Bagoly et al. 2009). As the real redshift distribution of 
GRBs cannot be changed by that much for the past few years, the only reason for it would 
be changing selection effects in the spectroscopic observations. 

These systematic biases will greatly affect various analyses. As long-duration GRBs 
come from the collapse of very massive fast-rotating stars with very short main-sequence 
lifetimes (Woosley & Bloom 2006), their rate density will provide a measurement of 
the massive star formation rate of our Universe. If we only deal with bursts with 
known spectroscopic redshifts, the derived star formation rate may be biased by the 
spectroscopic-redshift selection effects. The same problem comes with demographic studies 
of the GRB luminosity function. If only bursts with spectroscopic redshifts are used in 
constructing the luminosity functions, then the selection biases may distort the derived 
luminosity function. A third important problem is the recognition of the highest redshift 
bursts (say, with z > 7), against which the current spectroscopic methods may be heavily 
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biased. Swift is expected to have 5 — 10% of its bursts with z > 7 (Bromm & Loeb 2006), 
but the spectroscopic methods have only identified one such burst GRB090423, and that 
just a few month ago (Tanvir et al. 2009). 

A solution to these problems is to use the luminosity relations that work on luminosity 
indicators measured directly from the prompt gamma-ray emission. These luminosity 
relations are equations that connect a measurable property of the burst itself (the luminosity 
indicators) to the bursts energy (E-y or E^^iso) or the peak luminosity (L). The observed 
fiuence or peak flux can then be used (with the inverse-square law) to derive a distance 
to the burst and (for some fiducial cosmology) the burst redshift. The advantage of this 
procedure is that it apphes to all long class GRBs, not just to the minority selected to have 
a measured spectroscopic redshift. The disadvantage is that the uncertainties on the derived 
redshifts are much larger than those of the spectroscopic redshifts, and so the method can 
be used primarily for statistical or demographic purposes. The situation is similar to the 
case where photometric redshifts of distant galaxies in the Hubble Deep Field do not have 
the accuracy of the spectroscopic redshifts, yet nevertheless these photometric redshifts can 
be done for all the galaxies and are the cornerstone of all statistical studies of the field. 

The six luminosity relations we use in this article are: Tiag (spectral lag) - L relation 
(Norris, Marani, & Bonnell 2000), V (Variability) - L relation (Fenimore & Ramirez- Ruiz 
2000), Epeak (peak energy of the spectrum) - L relation (Schaefer 2003b), Ep^ak - Ej relation 
(also called Ghirlanda relation, Ghirlanda, Ghisellini, & Lazzati (2004)), trt (minimum 
rise time) - L relation (Schaefer 2002), and Np^ak (number of peaks in the light curve) - 
L relation (Schaefer 2002). It was claimed by Butler et al. (2007) that the Ep^ak - E^,iso 
relation (also called Amati relation, Amati et al. (2002)) differs somewhat between Swift 
and pre-Swift data, and hence the relation might be caused by the detection threshold 
effect of the instrument instead of the GRB itself. We tested this claim on the relations we 
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are using, as shown in Section 3. 

In this article, we will demonstrate a method to use these luminosity relations to 
estimate the GRBs redshifts. Here, we will only consider bursts with known spectroscopic 
redshifts, but our derivation of Zbest is based on the burst properties, and the effect of the 
known spectroscopic redshifts are negligible. This effect is tested later in the article. By 
applying our method to the bursts with accurate spectroscopic redshifts measured, we 
can compare our estimated redshifts with the measured spectroscopic redshifts. This will 
be the key test of our methods and the accuracy of our derived redshifts, as well as the 
applicability of our methods to bursts without spectroscopic redshift. With confidence in 
the reliability of the method, we can then apply it to all Swift bursts including those with 
no spectroscopic redshifts. Also, we will be able to apply this method to future bursts, and 
provide the community with a rapid notification about their redshifts. (Xiao & Schaefer 
2009) 



2. Improving the Luminosity Indicators 

In this paper, we adopted five luminosity indicators with six luminosity relations (with 
two relations for indicator Epeak- Epeak — L relation and Epeak — E^^iso relation). The 
details of the calculations for each of these indicators are described in Schaefer (2007). 
In this paper, we constructed the calculation independently, and made some significant 
improvements in the calculation of Tiag, trt and Np^ak- The five indicators are listed as 
below: 

(l)The spectral lag, Tiag, is the delay time between the soft and hard fight curves of 
a burst. By convention, we use the soft and hard energy bands to be those of BATSE 
channels 1 and 3, or Swift channels 2 and 4, covering roughly 25-50 keV and 100-350 keV. 
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For Tiag, by shifting the hard and soft hght curves of a GRB, and calculating the cross 
correlation between them, we are able to get a cross-correlation versus offset plot. The offset 
corresponding with the peak value of the cross-correlation is the lag time we need. The 
calculation is simple and easy for bright bursts, while for those faint ones, since the plot has 
significant scatters, the offset with the peak correlation (i.e. Tiag) is hard to evaluate under 
these noisy conditions. To find the offset when the cross-correlation achieves its peak value, 
we need to make a reasonable fit to the peak region of the cross correlation. As the shape 
and scatter of the plot varies from burst to burst, we cannot simply fit it with some specific 
function. If we fit it with a parabola, then any asymmetry in the cross correlation (as is 
often seen) will incorrectly shift the peak in the model by an amount depending on the 
range of offset included in the fit. And if we fit it with a high-order polynomial function, 
then the high-order terms will be unstable as they are trying to follow noise and regions far 
from the peak. What we did was to fit the cross correlation by polynomials with different 
orders (normally from 3 to 9), and choose the one which fits best in the very central region 
(around the peak) of the curve. Two examples are shown in Figure 1. 

A bootstrap procedure is previously used to calculate the uncertainties on the Tiag 
(Norris 2002). In our work, we are using a simple propagation method. The uncertainties on 
the cross-correlation amphtude points are calculated by simply evaluating the RMS scatters 
of these data points around the fitted curves. We are then able to generate the uncertainties 
for each of the fitting parameters, and the coefficient errors. Then the uncertainty for our 
Tiag value, ari^g , is generated by propagation of the uncertainties on the fitting parameters 
and the coefficients. 

(2) Variability (V) measures whether a light curve is spiky or smooth, and can be 
obtained by calculating the normalized variance of the original light curve around the 
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smoothed light curve. The calculation of V is as shown in Schaefer (2007): 

V =< [(C - Csmoothf - Oq]/ Csmooth,max > (l) 

where C is the count per time bin in the background subtracted light curve, with an 
uncertainty of ac- The time duration of the burst T„p is calculated as the summation of 
the time with the light curve brighter than 10% of its peak flux, and Cgmooth is the count in 
the smoothed light curve, with a box-smoothing width to be 30% of T^p. Csmooth,max is the 
peak value of C smooth- 

Prom equation (1), the 1 — a uncertainty of V can be propagated from the observational 
uncertainty of C in each time bin. To get rid of the cross-correlation effect between C and 
C smooth, for Q^ch Ci, we calculate it with 

^ i-l i+N/2 

Xi = Ci — Ci^smooth ~ ~ J^^ ^ ] Cj) (2) 

j=i-N/2 j=i+l 

where N is the box smooth bin (i.e. 30% of Tup). 

If wc neglect the uncertainty of ac, the 1 — a uncertainty of V can be propagated from 
equation (1) and the uncertainties of each Xi and Csmooth,max- 

The variability values (V) in Table 1 are about 10 times larger than those in Schaefer 
(2007). In this article, we are strictly following the definition and calculation of Equation 
(1). By checking the details of the calculation, we found that in Schaefer (2007), Cmax was 
used in the denominator instead of Csmooth,max, which caused the difference. The values 
of variability are very sensitive to the slightly change on its calculation function, and this 
might be one of the reasons for the large scatter of the variability-luminosity relation. 

(3) Epeak is the photon energy at which the uF^ spectrum is the brightest. By fitting 
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the GRB spectrum with a smoothly broken power law (Band et al. 1993) as 




(3) 



the peak energy Ep^ak and corresponding parameters a (the asymptotic power law index 
for photon energies below the break) and {3 (the power law index for photon energies above 
the break) can be obtained. Here $ is the usual differential photon spectrum {dN/ dE) as a 
function of the photon energy (E). A GRB spectrum is also able to be fit with a power-law 
with exponential cutoff model: 



with Epf.ak and power law index a being the fitting parameters. Here A and B are 
normalizing constants to indicate the brightness and constructed to ensure the continuity 
of the model spectrum. 

In this article, we are using E'peaifc values of the bursts from various of sources. All the 
values we used and the sources are hsted in Table 2. 

(4) The minimum rise time in the light curve, trt, was proposed for use in a luminosity 
relation by Schaefer (2002). The minimum rise time of a burst is taken to be the shortest 
time over which the light curve rises by half the peak flux of the pulse. In practice, especially 
for faint bursts with large Poisson noise, the rate difference between two close bins might 
be larger than half of the peak flux. As a result, we have to smooth the light curve before 
we calculate the rise time. The problem is then that if we smooth it too little, the apparent 
fastest rise time might be dominated by the Poisson noise, resulting in a too-small rise time, 
and if we smooth it too much, the smoothing effect will dominate, resulting in a rise time 
near the smoothing time bin. 

As the light curves vary greatly among different bursts, there is no specific box 
smoothing width that can satisfy a majority of bursts. Instead, we vary the box smoothing 



dN/dE = A £;"e-(2+°)^/-^p-'=, 



(4) 
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width from bins (that is, no binning) up to a relatively large number, say 50 bins, and for 
each of these smoothing widths, we generate a smoothed lightcurve from which a minimum 
rise time can be calculated. Of course, some of these light curves are over-smoothed and 
some are under-smoothed. In a minimum rise time versus smoothing width plot, we will 
have a monotonically rising curve, as shown in Figure 2. Although the shape of the curve 
varies amongst bursts, for most of the bursts there will be a region where the curve appears 
flat, or with a slightly increasing slope, which we call a plateau. It is easy to explain the 
existence of the plateau: in this region, the smoothing is enough that Poisson noise is 
negligible in determining the minimum rise times, while the smoothing is not so much that 
it determines the minimum rise time. On a plot of minimum rise time versus smoothing 
width, we can identify three regions: a fast rising region where the Poisson noise dominates, 
a nearly flat plateau region where we are seeing the real minimum rise time, following by 
another rising region for large box smoothing widths where the smoothing is dominating. 
We can make use of this plateau region by extrapolating it back to the zero-smoothing case 
(where the box smoothing width is zero), as shown in Figure 2. The intercept on the y-axis 
corresponds with the minimum rise time where the smoothing is effectively zero. With 
the extrapolation, we are sidestepping the regime where the Poisson noise dominates. As 
a result, the value of the intercept is just the minimum rise time we need, not affected by 
either the Poisson noise effect or the smoothing effect. 

For some of the extremely faint bursts, the Poisson noise dominant region and the 
smoothing effect dominant region will overlap with each other, and we are unable to flnd 
a plateau in the minimum rise time versus smooth width plot, for which the extrapolation 
cannot be made. Thus, our technique does not produce trt values for the faintest bursts. 

The uncertainty of the minimum rise time is calculated by simple propagation from 
the fltting parameters and the uncertainties on each individual point on the minimum rise 
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time versus smoothing width plot. The uncertainties on each of the individual points are 
dominated by the noise on each original data points in the hght curve. The noise on the 
peak flux will affect our criteria, by some factor of (Tc;/max(C), where C is the rate, and 
(7c is the uncertainty of the rate. In addition to that, random noise on the start and stop 
data point for each possible rise time will affect our determination, i.e. the real rise time 
between two data points may be larger /smaller than half of the peak flux, however, with 
the random noise on the start and stop points, we take it as our rise time, which is equal 
to half of the peak flux. This effect is also reflected on our criteria also, by a factor of 
~ 2 * (Tc'/max(C). As a result, for determine the uncertainties on each of the individual 
rise time on RT-smoothing width plot, we can change our criteria from 0.5 by a factor of 
2.25 * (Tc7/max(C), and record how much the resulted minimum rise time values changes. 
The uncertainties on the fitting parameters and the minimum rise time value can then be 
calculated from propagation. 

{5)Npgak is defined as the number of peaks in the light curve. With Cmax as the 
overall maximum of the background-subtracted light curve, we define a peak to be a local 
maximum that rises higher than Cmax/ 4: and is also separated from all other peaks by a 
local minimum that is at least Cmax/ 4: below the lower peak. In principle, Npeak is easy 
to count, either automatically or manually. In practice, we have the same problem of the 
Poisson noise and the smoothing factor effect, as what we had in the calculation of Tijr ■ 
Faint bursts will have their unsmoothed light curves dominated by apparent peaks produced 
by Poisson noise, resulting in large numbers of false peaks. A random noise spike can satisfy 
our definition for a peak if we dont smooth the light curve, yet if we smooth it too much 
there will always be just one peak. Here we adopted the same procedure of calculation as 
that in the calculation of trt ■ We vary the box smooth width from bins to a relatively 
large number, and calculate the number of peaks for each of the smoothed light curves. As 
in the case for trt , we see a fast falling curve (where Poisson noise is contributing spurious 
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peaks), with a plateau, where neither the Poisson noise effect nor the smoothing effect 
dominates. By extrapolating the plateau back to the y-axis, we get an Np^ak value for an 
unsmoothed case of the light curve, with the effects of Poisson noise removed. 

Each of the luminosity indicators discussed above has one or more corresponding 
luminosity relations. These relations are Tiag — L, V — L, Ep^ak — L, Ep^ak — E^ (so called 
Amati relation) and Epeak — E^^iso (so called Ghirlanda relation), trt — L, and Npeak — L 
relation, which are described and explained in details in Schaefer (2007). In our redshift 
calculation program, we are not including Amati relation for the following two reasons: 
First, the Amati relation has been challenged as it returns ambiguous redshifts (Li 2007), 
while the Ep^ak — L relation (also the Tiag — L, V — L, trt — L, and Np^ak — L relations) 
passed the same test (Schaefer & CoUazzi 2007). Second, the physics of the Amati relation 
is nearly the same as that of Ghirlandas Epeak — E^ relation, except that Ghirlandas 
Epeak ~ relation has a correction for the jet opening angle, making it much tighter. 
There is another luminosity relation which relates Epeak, an effective duration (called T45), 
and the luminosity by Firmani et al. (2006). However, with more GRBs and Epeak data 
added in, it is realized that this luminosity relation is not making any improvement on the 
Epeak — L relation (CoUazzi & Schaefer 2007). 

3. Data for 107 Long GRBs and Six Luminosity Relations 

We have 107 long GRBs with their spectroscopic/photometric redshifts measured, 
ranging from Feb. 28, 1997 (GRB970228) to July 21, 2008 (GRB080721), observed by 
BATSE, Konus, HETE, and Swift. The average redshift for pre-Swift bursts is about 1.50 
and that of the Swift bursts is 2.15. 

The light curves of Swift GRBs are generated from the original data published 
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on the legacy ftp sitejl], and the Swift Software ver 2.9 (HEAsoft 6.5). To gener- 
ate a background-subtracted hght curve of a GRB, we downloaded an event file 
swOOxxxxxxOOO.bevshsp_uf.evt.gz and a mask file swOOxxxxxxOOObcbdq.hk.gz, with xxxxxx 
be the six digit Swift trigger number. By running a task 'batbinevt', we can specify the 
time interval, energy bins, time bin method, output file name and format on generating the 
light curve. In our work, we adopted a time interval of 0.064 s with uniform time bins, and 
four continuous energy bands (15-25 keV, 25-50 keV, 50-100 keV and 100-350 keV). For the 
calculation of V, t/jt and Np^aki we have been using the light curve over the whole energy 
range (15-350 keV), and for the calculation of Tiag value, we use 25-50 keV and 100-350 keV. 
The light curves of pre- Swift bursts are obtained from our previous work (Schaefer 2007). 

We calculated the Tiag, V, trt, and Npeak values for each of these bursts, as listed 
in Table 1. The first column lists the ID number of GRBs. The second column lists the 
satellite with the detection of the burst. Column three to six show all the calculated 
indicator values, with the name of the indicators shown on the header row. Indicators not 
measured due to low signal-to-noise ratio of the burst are represented as 

The values of Ep^ak as well as the power law indexes a, (3 in Band's smoothly broken 
power law model (Band et al. 1993) or Epeak and a values in a power law with exponential 
cutoff model are obtained from various sources, as shown in Table 2. The first column in 
Table 2 lists the ID number of the GRBs. The second column lists the satellite with the 
detection of the burst. Column three to six are the values of -Epeafc and the power law index 
a and (3 values, as well as the reference sources of these values. Our jet break time {tjet) 
values from optical observations and their sources are also listed in Table 2, column seven 
and eight. Bursts without measured jet break time are filled by '. . .' also. From the table 
we see that only 33 of the bursts have their optical tjet reported. Various jet break time in 
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the X-ray afterglows have been reported, however, as there are usually multiple breaks for 
the X-ray afterglows, which are not well understood and distinguished for their causes, we 
are not including any of these reported tjet values from X-ray detections. Values in square 
brackets for a and (3 are assumed values for those bursts without exact a and (3 values 
measured, which is taken to be the average value of known a and j3 (Schaefer et al. 1994; 
Krimm, et al. 2009; Kaneko et al. 2006; Band, et al. 1993). Some uncertainties of -Epeafc 
and tjet are also quoted in square brackets. These uncertainties are assumed conservative 
assumed values, which are normally 10% of the measured Ep^ak or tjet values. All of the 
peak flux (P) and fiuence (S) values of pre-Swift bursts we use here are as same as what 
were used in Schaefer (2007), and those of Swift bursts are from the data table on the Swift 
webpagcl. All the quoted error bars in Table 1 and Table 2 are converted to 1 — cr level. 

The luminosity relations are all expressed as power laws, and we can make a linear 
fit on the logarithms of the redshift-corrected luminosity indicators and the logarithms of 
the burst luminosities. In Figure 3, we display the data and the best fits for each of these 
luminosity relations, for the concordance cosmological model (w = -1, and Qm = 0.27 in 
a fiat universe). As there are significant uncertainties in both the horizontal and vertical 
axes, and some intrinsic scatters independent of both the luminosity and the redshift, we 
performed the ordinary least squares without any weighting (Isobe et al. 1990, Schaefer 
2007). As both the indicator and luminosity are caused by some other parameter (like the 
jet Lorentz factor in particular), and so both indicator and luminosity are correlated. In 
this case, we assume that these two variables in the luminosity relations are not directly 
causative, and the bisector of the two ordinary least squares (Isobe et al. 1990) has been 
used. More details of the fitting process are referred to Schaefer (2007). As Ep^ak plays an 
important role in the calculation of Pboio, the uncertainties on luminosity is then correlated 



^http: / / swift.gsfc.nasa.gov/docs/swift / archive/grb_table/ 
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with the uncertainties on Epeak- In our work, we ignored that effect, and simply assumed 
that Phoio and Ep^ak have uncorrelated errors. The best fitting function for each of the 
relations are shown in Table 3. 

The V — L relation is rather scattered, as we can see from Figure 3 and Table 3. The 
scatter is larger than one could expected for any linear relations. As the calculation of 
V is not well defined, and V — L relation is our most noisy luminosity relations. In this 
case, we decided not to include V — L relation in our later calculation for the redshifts, 
although we listed all the calculated V values for the bursts in our Tables. As a result, the 
maximum number of luminosity indicators that we can use is 4, and the maximum number 
of luminosity relations is 5. 

Butler et al. (2007) estimated the Amati's relation {Ep^ak ~ E^,iso) with Swift and 
pre-Swift data. They report an inconsistency in the relations from the two data sets, which 
they then attributed to differences in threshold between Swift and earlier detectors. Their 
claimed difference has not been reproduced by other groups (including Cabrera et al. 2007; 
Schaefer 2007b; Krimm et al. 2009), while their claimed threshold effects have been found 
to not significantly affect the observed Amati relation (Schaefer 2007b; Nava et al. 2008; 
Ghirlanda et al. 2008; but see Shahmoradi & Nemiroff 2009). Indeed, their analysis is 
based on Bayesian priors which systematically push high-£^peafe values below ~ 400 keV, 
as demonstrated by detailed comparisons with Konus, Suzaku, and RHESSI measures. 
Nevertheless, in this paper we can perform yet another test to see whether the claimed 
threshold differences between Swift and pre-Swift bursts cause any significant change in 
the luminosity relations. Butler et al. (2007) paper is including Swift BAT bursts between 
GRBs 041220 and 070509, 77 of which have spectroscopic redshift measured. While in our 
analysis on Swift GRBs, we are including all bursts with spectroscopic redshift between 
050126 and 080721. With both samples being based on largely overlapping samples selected 



-16- 



in a nearly identical manner, we conclude that the flux limits of the two samples are 
essentially identical. 

To this end, we have separately fitted the pre-Swift and Swift data. This test was 
only done for four luminosity relations {riag — L, V — L, Epeak — L, and trt — L), with the 
Epeak — relation having too few bursts, and the N^^ak — L relation not being usable for 
the comparison as a limit. The best fit luminosity relations are given in Table 4. The data 
and best fit models are displayed in Figure 4. At first glance, we see that the difference 
between the two best fit lines are small compared to the scatter in the data, and a detailed 
analysis is described below. 

We made a F-test for the fitting results for these four relations. First we made a 

bisector linear fit on the combined data with both pre-Swift and Swift bursts, and recorded 

the value of the fit as Xjoint- Then we separate the data to two sample sets, pre-Swift 

and Swift, and made the same bisector linear fit separately on each set of the data. The 

sum of the two values for the separately fitted lines are recorded as Xseparate- Then 

the F value can be calculated as 

jp _ Xjoint/ {^pre + ^sw ~ 2) 
Xseparate/ i^pre + ^sw — 4) 
where Np^e is the number of pre-Swift bursts, and Ngw is the number of Swift bursts. 

Npre + Ns^ — 2 is the degree of freedom of the fitting on the combined data, and Npj.e + Ng^ — 4 

is the degrees of freedom of the fitting on separated data. 

These F values for each of the luminosity relations are listed in Table 4. If the pre-Swift 
and Swift relations differ much from each other, the separate fitting have been significantly 
improved over the fitting on the mixed data, the F value would be much larger than unity. 
Otherwise, if there is no significant difference between pre-Swift and Swift relations, the F 
value would be around unity. From Table 4 we see that F is rather close to unity. All of the 
results shows that the separately fitted result is not significantly improved over the fitted 
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results on all data mixed together, which means that luminosity relations for pre-Swift 
luminosity relations and Swift luminosity relations do not significantly differ from each 
other. And by looking at the plots in Figure 4, we see that (1) the envelope of squares and 
diamonds are indistinguishable and (2) the pre-Swift and Swift best fit lines are close to 
each other compared to the scatter in the data. And hence, we have no significant evidence 
that these four luminosity relations differ for Swift bursts. 

The 1 — a range on the normalization difference between Swift and pre-Swift bursts for 
all these four luminosity relations are also listed in Table 4. We can make an analysis with 
the normalization difference on Amati's relation claimed by Butler et al. (2007), which is 
corresponding with a 0.39 difference in log space. By making the comparison between the 
Butler's factor (0.39) and our normahzation difference, we can exclude Butlers factor at a 
2.5 sigma level for Tiag — L relation, a 3.7 sigma level for trt — L relation, a 1.3 sigma level 
for Epeak ~ L relation (which is not significant), and we cannot exclude the Butlers factor 
for y — L relation. 

4. Method for Calculating Redshift 

Our aim in this paper is to test our method of redshift calculation. We applied it to 
the bursts with known spectroscopic redshifts {zspecji and if it works well, we will be able to 
apply it to all the long GRBs in our future work. Although we are dealing with the bursts 
with known Zgpec-i these Zgpec are only involved in the fitting of luminosity relations, and this 
effect is negligible in our calculation. Only after we have derived our redshift Zbest from the 
luminosity relations will we compare them with the known Zgpec to test our accuracy. 

Below we will describe how our method applies on one Gamma-Ray Burst step by step: 

(1) First, we measure each of the luminosity indicators of the burst. The definition and 
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method of calculations have been discussed in Sections 2 & 3. The results of the indicators 
for all bursts in the sample are hsted in Table 1 & 2. 

(2) We next derive the luminosity values for each relations from Table 3. A complexity 
is that the luminosity relations depend on the redshift of the burst (so as to correct the 
luminosity indicators back to the burst rest frame), so we have to perform this calculation 
for an array of trial redshifts (we take it to be from redshifts of to 20 at intervals of 0.005), 
and then we will obtain a list of luminosities (or isotropic energies for Ghirlanda's relation) 
values depending on the list of trial redshifts for each of the indicators. We notate each of 
these calculated luminosities (as a function of redshift Zfriai) for the ith relation as Li{ztriai) 
(or E^^i for the Ghirlanda's relation) . 

(3) With the values of the peak flux P, fluence S, Ep^^k and the power law indexes in 
the broken power law model a and (5 (or a from the power law with exponential cutoff 
model), the bolometric peak flux Phoio and fluence Sboio can be calculated. The range for 
'bolometric' is set to be 1 keV to 10000 keV in the ORB rest frame, and the equations for 
the detailed calculation are referred to Schaefer (2007). As a result, for each burst, we 
calculated Pboioi^triai) and Sboioi^triai) for each trial redshift value from to 20. 

For those bursts with tjet values, the jet opening angle 9jet (in units of degrees) can be 
calculated as 

e,et = 0.161 (6) 

where tj^t is the jet break time in the unit of days, n is the density of the circumburst 
medium in particles per cubic centimeter, rj^ is the radiative efficiency, and E^^iso,52 is the 
isotropic energy in units of 10^^ erg (Rhoads 1997 & Sari et al. 1999). We simply adopt 
77-y = 0.2 and n — 3cm~^ in equation 6. The beaming factor Fheam, is then calculated as 

Fbeam ^ cos e jet- (7) 

Prom above we see that, since both Oj^t and E^^iso are redshift sensitive, our calculated 
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Fbeam value also varies between different Ztriai values. 

(4) From all the parameters above, for each of the indicators, a list of the luminosity 
distances can be calculated as 



For Ghirlanda's relation, the hst of luminosity distances is calculated as 



7 / s E^^i(^Ztrial^{^ ~l~ Atrial] /f^\ 

aL,i[^trial) = \ -r-^ ^ 7- ^- (9; 

y ^'Krbeam^boloyZtrial) 

Then from each list of luminosity distance above, the distance modulus can be obtained: 

l^i{Ztricd) = 5 l0g[dL,i{ztrial)] " 5, (10) 

with dL^i{ztriai) expressed in units of parsecs. The uncertainties are propagated strictly 
following the calculation. 

From all above, we are able to get up to four lists of measured distance moduli iJ>i{ztriai) 
with their 1 - a uncertainties: Hn,^ ± a^,^^^, fiE,,^>, ± ct^e^^^,, ^E^.ak-E, ± ^he^^^^-e,^ and 
lirET '^t^Tf^j.- As have asymmetric uncertainties for Epeak, we carry the uncertainties 
on both directions in the calculation, and generated both plus and minus uncertainties for 
I^Ep^ak I^Epeak-E-y- I^Npeak ^Iso calculatcd, as a lower limit on the distance modulus. All 
of these distance moduli are a function of the assumed Ztriai for < Ztriai < 20. 

(5) Given each trial redshift, we can calculate its distance modulus directly from the 
cosmo logical model, iicos{ztriai) ■ Here we adopt the concordance model, with equation 
of state for dark energy p = wpc^, w = —1, VIm = 0.27, fl\ = 1 — VIm = 0.73, and 
Hq — 70{km/s)/Mpc. In this case, the luminosity distance can be expressed as 

dLiZtrial) = CH^\1 + Ztriai) f '"^^ dz'[{l + z'f^lu + ^kY^''' ■ (H) 
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Prom the equation above and our list of trial redshifts, a list of luminosity distances 
dL{.Ztriai) will be calculated, and also a list of distance modulus fJ^cosi^triai) which equals 
5 log[dL{ztriai)] — 5. The Hcosi^triai) valucs are only depending on the trial redshift (running 
from to 20) and the cosmological model we choose. 

(6) For each of the trial redshifts from to 20, we have distance moduli lists of A*T(„g) 
A'Speofc) A*Epeafc-£^7' I^trt ^^^ug wlth thclr 1 — a uncertainties as well as /icos- We can then 
compare these iXi{ztriai) with Hcosiztriai) in a sense. Thus, „ = [(A*r,,, - IJ'cos) / t^r, f and 
so on for the other relations. We get a versus Ztriai plot, as shown in the left panel of Figure 
5. Then sum over all the xl we get a xltai ixLai = xI,^^ + X%^^, + X%^^,-e, + xI^^ + X%^^J 
versus Ztriai plot , as in the right panel of Figure 5. Our best redshift (zbest) corresponds 
with the minimum Xtotah where the luminosity relations and the cosmological model agree 
with each other best. We are also able to find the uncertainties of our Zf,est- By searching 
through the xLtai - ^triai plot, we can find the redshifts with which the xLtai = xlotai,min + 
which corresponds with the edges of the 1 — a range of our redshift. Similarly, the redshifts 
with Xtotai = Xtotai,min + ^ gives US the 2 - cT range and that with Xtotai = Xtotai,min + 9 gives 
the 3 — 0" range of our redshift. 

5. Results 

The resulting values of Zb^st and 1 — a range of z for each burst are shown in Table 5. 
We have also collected the spectroscopic redshifts for these bursts (see values and references 
in Table 5). The last column of Table 5 lists the effective luminosity relations wc used for 
each burst in our calculation. The number of luminosity indicators used in the redshift 
calculation for each of the burst is fisted in Table 5. A comparison can be made between 
our calculated redshifts and their spectroscopic (or photometric) redshifts. The comparison 
plot is shown is Figure 6, and below are some of the conclusions after we made the analysis: 
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(1) Of the total 115 bursts, 8 have only a lower z hmit that can be calculated. For 
the remaining 107 bursts, we took the spectroscopic redshifts Zspec as the model value, and 
our calculated redshift Zhest as a measured value, with the 1 — a uncertainty of cr^^^^^ and 
(72,.. . Then the can be calculated as 



with aj equals to CFj^zpi^s ^j,ZMinusi depending on whether our Zi^^est is smaller or larger 
than the Zi^spec- Each j represents an index number identifying the burst. The number of 
degrees of freedom in this comparison equals the number of bursts [Ni^^st) , so the reduced 
is X^/^burst- The reduced value is 1.28, which is somewhat larger than unity. While 
after excluding one 3 — cr outliers GRB010222, whose contribution to is as high as 44, our 
reduced is equal to 0.86. This is certainly not a significant deviation from unity. So we 
conclude that the scatter in Figure 6 is consistent with our quoted error bars being correct. 

(2) Of the 107 bursts with their Zf,est calculated, 73 has their Zgpec falling into the 1 — cr 
range of our z. The ratio of the numbers is about 70%, which is slightly larger than the 
ideal case 68.5%. And for the 8 bursts with only lower z hmit calculated, 6 of them have 
their Zgpec larger than our 1 — a lower limit of z. This is another way of testing our quoted 
error bars, and as in the previous item, we find no significant deviation from the expected 
results. As such, to a close degree, we see that our derived error bars are accurate. 

(3) We can test to see whether our z^est is biased high or low. For this, we calculated 
the average value of logio{zbest/ ^spec)- If our result is unbiased, the average should be zero 
to within the error bars. We find the average is 0.01. This demonstrates that our Zbest is 
not biased to within the 1% level. 

(4) To test the accuracy of our Zbest comparing with the Zgpec: we calculated the RMS 
scatter of logio{zbest/ ^spec)- The result comes out to be 0.26, while the RMS scatter of 




(12) 
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logio{zspec) is 0.30. In Schaefer (2007) it was pointed out that the accuracy of the redshift 
estimation is 26% (corresponding to a logio{z{,est/zspec) RMS of 0.11), which is better than 
what we are claiming here. The reason for the larger RMS scatter is, although we are 
dealing with the GRBs with known Zspec, we are not making any use of the Zgpec in our 
whole calculation, and our luminosity L, luminosity distance dL, distance modulus Hind, 
as well as the bolometric flux and fiucncc Pboio and Sboio are all varying with our trial 
redshift. This brought us one extra degree of freedom in the calculation, which caused 
larger uncertainties in our result. When the known Zspec value is used (as for the Hubble 
diagram work in Schaefer 2007), the scatter becomes substantially small as compared to our 
work in this paper. Another reason for the large RMS is because it is mostly dominated by 
the very noise bursts (those bursts with low signal to noise ratio, inaccurate measurement 
of luminosity indicators, and large error bars for z^est results), which are not able to be used 
for the notification of the redshifts. If we make a selection on bursts with relatively accurate 
redshift estimation, say bursts with cr, . /zu * < 0.5 and cr, , , < 1, we get the RMS of 
log{zi,est/ Zspec) of 0.19, which is much smaller, and the RMS scatter of logiQ{zspec) for this 
subsample is 0.28. 

(5)We made the same calculation with pre-Swift luminosity relations on calculating 
Swift redshifts and Swift luminosity relations on calculating pre-Swift redshifts, for which 
we call it Zhest,ii- Our calculation on Zhest,ii is totally independent on the spectroscopic 
redshifts, as for each burst, the luminosity relations used in the calculation are calibrated 
independent of the Zspec for any individual burst. The comparison between Zbest and 
Zhestji will also show the difference between Swift and pie-Swift luminosity relations. The 
comparison plot between Zspec and Zbestji is shown on figure 7. From the comparison 
between figure 6 and figure 7 we see that, the scatter and distribution of Zbest and Zbestji do 
not differ significantly from each other. Actually, from our calculation, the average value 
of logio{zbest,ii/zspec) is -0.02, and the RMS scatter of logio{zbest,ii/ Zspec) is 0.27, both of 
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which are equal to those of our Zbest value within error bars. This result demonstrate that 
those two sets of Zbest values do not differ from each other, which means that these two sets 
of luminosity relations do not have significant difference. It also tells us that the effect of 
redshift involved in our calculation (in fitting luminosity relations) are negligible. 

(6)We need to verify whether our result is effective in selecting high redshift bursts. 
If our predicted redshift is z, the possibility of a real redshift to be higher than z and 
lower than z are both 50%, which cannot be used as a test. However, by considering our 
uncertainties of Zbest, if our predicted redshift is 2z, we can make a test by counting how 
many of the GRBs are with real redshifts larger than z. As there are not many bursts with 
high redshifts, a test is done on a relatively lower redshift region, where most of the GRBs 
are involved. We picked up all our GRBs with predicted z > A, the total number is 12, and 
out of 12 have their spectroscopic redshift z < 2. This result tells us that our method is 
actually effective in demographic studies and in picking up high z bursts, if we take into the 
consideration of the error bars before we make the prediction. 

Prom all the analysis above, we can conclude that our Zbest is not biased on average, 
and our 1 — a error bars are accurate. We can claim that our method works well on the 
bursts with known spectroscopic redshifts, and can be applied to all long GRBs (even 
without their Zgpec measured). 

6. Conclusions 

In this paper, we developed a method to calculate the redshifts for long GRBs, using 
their light curves and spectra. We applied our method to bursts with known spectroscopic 
redshifts, detected by BATSE, HETE, Konus and Swift. By comparing our calculated 
redshifts with their spectroscopic redshifts, we are able to examine the accuracy of our 
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method. 

We compared each of the luminosity relations for pre- Swift and Swift bursts by making 
a F-test. With the F values close to unity, we have significant evidence against any claim 
that the relations are caused by the detection threshold effects or any other artificial effects 
of the instruments. 

We compared our results with the spectroscopic redshifts. We find that our Zh^st are 
not biased (with the average value of logiQ^z^estl ^spec) equal to 0.01), and our reported 
1 — (7 error bars are good (with x^ed — 0-86, and 70% of the Zgpec fall into the 1 — a region 
of Zbest)- Our accuracy on the redshifts are not as accurate as those from spectroscopy, 
yet nevertheless with a reasonable accuracy for demographical and statistical studies, with 
the RMS of logio{zbest/ ^spec) is 0.26. The RMS value is about twice as what was found 
in Schaefer (2007). One of the reason is, in Schaefer (2007), the accuracy is calculated 
assuming a known Zspec, and in this paper, as we are treating the unknown-redshift case, so 
extra degrees of freedom has been brought in the calculation, which caused the accuracy to 
get worse by about a factor of 2. Another reason is that the large RMS is dominated by 
those faint and noisy burst, and for a subsample with a. /x, ^ < 0.5 and o", , , < 1, 
we get the RMS of [log^Zi^est/ Zgpec)] of 0.19, which is much smaller. As our ^r^est are from the 
light curves, the spectra and the concordance cosmological model, it is independent of the 
spectroscopic redshift. As a result, our method can be applied to all long GRBs. 

For Swift bursts, as we are measuring over a relatively narrow energy band (15 keV - 
350 keV), the uncertainties in the calculation of peak energy Ep^^k is large, and it becomes 
one of the main restrictions on our accuracy. With the launch of Fermi, we are hoping to 
get bursts with more accurately measured Ep^ak values and light curves covering a broader 
energy band. We are expecting a substantial improvement in the accuracy of redshifts for 
Fermi bursts. 
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For the next step, we will apply our method to all Swift long GRBs, aiming to get a 
nearly-complete Swift GRBs redshift catalog. Such a catalog will inevitably be incomplete 
due to bursts with incomplete light curves and bursts too faint for their properties to 
be usefully measured. Our resulting redshifts will have an accuracy worse than those 
obtainable with optical spectroscopy, yet our accuracy will be good enough for various 
important statistical studies. As such, our catalog will be used for the demographic studies, 
without the detection threshold effect of the spectroscopic redshift measurements. We are 
also trying to deal with all the future bursts, and to provide a rapid notification of the 
redshift on the GCN circular to the community. We are hoping to find some possible high 
redshift (e.g. z > 7) and possible low redshift (e.g. z < 0.3) GRBs, which will be important 
in the observation of the high redshift universe and the GRB-SNe connection study. 



This work is supported by NASA under grant NESSF 07Astro07F-0029. 
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Table 1. Luminosity Indicators. 



GRB Satellite Tiag{sec) V t^t (sec) Np^ak 



970228 


Konus 







016 


± 





010 












2 


970508 


BATSE 


0.49 ± 0.02 





018 


± 





004 





65 


± 





07 


1 


970828 


Konus 







052 


± 





005 





36 


± 





14 


4 


971214 


BATSE 


0.03 ± 0.05 





048 


± 





002 












2 


980703 


BATSE 


0.69 ± 0.02 





024 


± 





001 


3 


00 


± 





19 


1 


990123 


BATSE 


0.07 ± 0.01 





059 


± 





003 












3 


990506 


BATSE 


0.04 ± 0.01 





337 


± 





001 





13 


± 





01 


12 


990510 


BATSE 


0.03 ± 0.01 





118 


± 





001 





13 


± 





01 


8 


990705 


Konus 







097 


± 





004 





62 


± 





37 


4 


991208 


Konus 







023 


± 





003 





27 


± 





01 


4 


991216 


BATSE 


0.03 ± 0.01 





062 


± 





003 





09 


± 





01 


5 


000131 


Konus 







056 


± 





005 





84 


± 





39 


2 


000210 


Konus 







018 


± 





002 





45 


± 





03 


1 


000911 


Konus 







122 







013 





07 


± 





22 


5 


000926 


Konus 







326 


± 





034 












4 


010222 


Konus 







143 


± 





004 





45 


± 





01 


6 


010921 


HETE 


1.00 ± 0.04 





008 


± 





006 


4 


31 


± 





71 


1 


020124 


HETE 


0.07 ± 0.06 





266 


± 





040 





59 


± 





17 


3 


020405 


Konus 







104 


± 





007 





48 


± 





09 


2 


020813 


HETE 


0.15 ± 0.01 





164 


it 





004 





59 







05 


5 


021004 


HETE 


0.71 ± 0.19 





035 


± 





067 


1 


23 


± 





96 


2 


021211 


HETE 


0.31 ± 0.01 





006 


± 





003 





57 


± 





01 


1 


030115 


HETE 


0.44 ± 0.06 





020 


± 





020 





70 


± 





40 


1 


030226 


HETE 


0.31 ± 0.22 





033 


± 





029 


1 


76 


± 


1 


15 


3 


030323 


HETE 







021 


± 





338 












2 


030328 


HETE 


0.08 ± 0.08 





024 


± 





003 


1 


69 


± 





81 


2 


030329 


HETE 


0.15 ± 0.01 





065 


± 





002 





66 


± 





01 


2 


030429 


HETE 


0.03 ± 0.17 





220 


± 





135 












2 


030528 


HETE 


12.56 ± 0.14 





017 


± 





010 


2 


13 


± 





42 


1 


040924 


HETE 


0.90 ± 0.01 





060 


± 





003 





33 


± 





17 


1 


041006 


HETE 







050 


± 





002 


1 


28 


± 





01 


3 


050408 


HETE 


0.31 ± 0.02 





082 


± 





005 





49 


± 





02 


1 
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Table 1 — Continued 



GRB Satellite Tias(sec) V rjiT(sec) Np^ak 



051022 


Konus 




0.088 ± 0.008 





19 


± 





04 


1 


050126 


Swift 


2.74 ± 0.02 


-0.010 ± 0.065 


1 


58 


± 


1 


91 


1 


050223 


Swift 




0.111 ± 0.094 












1 


050315 


Swift 




0.032 ± 0.016 


1 


97 


± 


1 


62 


2 


050401 


Swift 


0.06 ± 0.02 


0.187 ± 0.019 





25 


± 





16 


3 


050406 


Swift 




0.020 ± 0.274 












2 


050416A 


Swift 




0.021 ± 0.030 





54 


± 





06 


1 


050505 


Swift 


0.71 ± 0.13 


0.076 ± 0.031 





60 


± 





21 


3 


050525A 


Swift 


0.12 ± 0.01 


0.093 ± 0.003 





32 


± 





01 


2 


050603 


Swift 


-0.01 ± 0.01 


0.125 ± 0.014 





19 


± 





01 


1 


050730 


Swift 




0.027 ± 0.066 












2 


050802 


Swift 




0.070 ± 0.036 


2 


03 


± 


1 


02 


4 


050814 


Swift 




-0.009 ± 0.180 












2 


050820A 


Swift 




0.061 ± 0.033 


1 


01 


± 





75 


3 


050824 


Swift 




0.289 ± 0.640 












1 


050826 


Swift 




0.063 ± 0.105 


1 


11 


± 


2 


28 


1 


050908 


Swift 




-0.017 ± 0.046 


1 


10 


± 


1 


47 


1 


050922C 


Swift 


0.06 ± 0.01 


0.015 ± 0.003 





13 


± 





01 


2 


051016B 


Swift 




0.008 ± 0.030 












2 


051109A 


Swift 




-0.006 ± 0.025 





70 


± 


1 


25 


1 


051111 


Swift 


1.70 ± 0.07 


0.009 ± 0.004 


1 


80 


± 





24 


1 


060108 


Swift 




0.006 ± 0.040 












2 


060115 


Swift 




0.019 ± 0.029 


1 


11 


± 


1 


71 


2 


060206 


Swift 


0.01 ± 0.03 


0.007 ± 0.004 


1 


16 


± 





18 


1 


060210 


Swift 


0.15 ± 0.17 


0.183 ± 0.033 





73 


± 





50 


4 


060223A 


Swift 




0.0.36 ± 0.021 





41 







23 


4 


060418 


Swift 


0.22 ± 0.03 


0.104 ± 0.008 





67 


± 





08 


2 


060502A 


Swift 


4.90 ± 0.11 


0.004 ± 0.010 


2 


94 


± 


1 


19 


1 


060510B 


Swift 




0.110 ± 0.060 












4 


060512 


Swift 




0.043 ± 0.173 












1 


060522 


Swift 




0.034 ± 0.185 












1 


060526 


Swift 


0.17 ± 0.09 


0.085 ± 0.030 





38 


± 





11 


2 
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Table 1 — Continued 



GRB Satellite Tias(sec) V rjiT(sec) Wpeafc 



060604 


Swift 




0.080 ± 0.338 








2 


060605 


Swift 




-0.013 ± 0.068 


1.22 


± 


0.72 


3 


060607A 


Swift 


1.98 ± 0.11 


0.025 d= 0.008 


1.23 


± 


0.68 


1 


060707 


Swift 




0.050 ± 0.054 








2 


060714 


Swift 




0.125 ± 0.022 








^3 


060729 


Swift 




0.092 ± 0.041 








2 


060814 


Swift 


0.29 ± 0.03 


0.040 ± 0.003 


1.65 


± 


0.24 


2 


060904B 


Swift 


0.36 ± 0.09 


0.003 ± 0.008 


1.00 


± 


0.16 


1 


060908 


Swift 


0.26 it 0.06 


0.061 ± 0.008 


0.52 


± 


0.09 


3 


060926 


Swift 


1.03 ± 0.11 


0.148 ± 0.050 








2 


060927 


Swift 


0.12 ± 0.04 


0.094 ± 0.010 


0.46 


± 


0.12 


2 


061007 


Swift 


0.11 ± 0.01 


0.066 ± 0.003 


0.38 


± 


0.02 


4 


061110A 


Swift 




-0.038 ± 0.050 








1 


061110B 


Swift 


0.24 ± 0.36 


0.155 ± 0.064 


0.79 


± 


0.64 


9 


061121 


Swift 


0.03 ± 0.01 


0.050 ± 0.003 


0.98 


± 


0.19 


2 


061222B 


Swift 




0.024 ± 0.043 








2 


070110 


Swift 




-0.010 ± 0.031 








1 


070208 


Swift 




0.083 ± 0.211 








2 


070318 


Swift 




0.037 ± 0.008 


0.72 


± 


0.24 


1 


070411 


Swift 




0.041 ± 0.029 








2 


070506 


Swift 


2.52 ± 0.04 


0.010 ± 0.030 


0.12 


± 


0.06 


1 


070508 


Swift 


0.04 ± 0.01 


0.106 ± 0.003 


0.20 


± 


0.01 


4 


070521 


Swift 


0.04 ± 0.01 


0.116 ± 0.004 


0.58 


± 


0.06 


5 


070529 


Swift 




0.170 ± 0.091 








1 


070611 


Swift 




0.053 ± 0.080 








1 


070612A 


Swift 




0.032 ± 0.023 


2.49 


± 


1.48 


2 


070714B 


Swift 


0.03 ± 0.01 


0.164 ± 0.021 


0.45 


± 


0.04 


1 


070802 


Swift 




-0.156 ± 0.150 








1 


070810A 


Swift 


1.09 ± 0.23 


-0.006 ± 0.015 


0.73 


± 


0.22 


1 


071003 


Swift 


0.38 ± 0.05 


0.072 ± 0.007 


0.88 


± 


0.07 


4 


071010A 


Swift 




-0.076 ± 0.153 








1 


071010B 


Swift 


0.84 ± 0.04 


0.010 ± 0.003 


1.21 


± 


0.03 


1 



Tabic 1 — Contimicd 



GRB Satellite Tiagisec) V t^t (sec) A'^peafe 



071031 Swift ... -0.038 ± 0.108 ... 2 



071117 


Swift 


0. 


,60 


± 


0, 


,01 


0.009 


± 


0.003 


0, 


,20 


± 


0, 


,02 


1 


071122 


Swift 












0.391 


± 


0.392 












1 


080210 


Swift 


0, 


,53 


± 


0, 


,17 


0.019 


± 


0.013 


0, 


,57 


± 


0, 


,44 


3 


080310 


Swift 












0.038 


± 


0.021 


0, 


,41 


± 


0, 


,55 


3 


080319B 


Swift 


0. 


,02 


± 


0, 


,01 


0.031 


± 


0.003 


0, 


,14 


± 


0, 


,01 


10 


080319C 


Swift 












0.042 


± 


0.007 


0, 


,21 


± 


0, 


,12 


4 


080330 


Swift 












0.109 


± 


0.060 












3 


080411 


Swift 


0, 


,21 


± 


0, 


,01 


0.167 


± 


0.003 


0, 


,65 


± 


0, 


,01 


2 


080413A 


Swift 


0, 


,13 


± 


0, 


,03 


0.078 


± 


0.004 


0, 


,23 


± 


0, 


,03 


3 


080413B 


Swift 


0. 


,23 




0, 


,01 


0.004 


± 


0.003 


0, 


,50 


± 


0, 


,03 


1 


080430 


Swift 


0. 


,68 


± 


0, 


,08 


0.009 


± 


0.004 


0, 


,76 


± 


0, 


,12 


1 


080516 


Swift 


0. 


,15 


± 


0, 


,01 


0.168 


± 


0.055 












2 


080520 


Swift 












0.037 


± 


0.098 












1 


080603B 


Swift 


0, 


,08 


± 


0, 


,01 


0.283 


± 


0.010 


0, 


,22 


± 


0, 


,03 


6 


080605 


Swift 


0, 


,11 


± 


0, 


,01 


0.057 


± 


0.003 


0, 


,22 


± 


0, 


,01 


4 


080607 


Swift 


0. 


,04 


± 


0, 


,01 


0.035 


± 


0.003 


0, 


,18 


± 


0, 


,06 


6 


080707 


Swift 












0.093 


± 


0.032 












2 


080721 


Swift 


0. 


,13 


± 


0, 


,05 


0.048 


± 


0.009 


0, 


,09 


± 


0, 


,04 


4 



Table 2. Luminosity Indicators. 



GRB 


Satellite 




a* 




Ref.t' 


tjetiday) 


Rel 


970228 


Konus 


11 5+38 


-1.54 + 0.08 


-2.5 ± 0.4 








970508 


BATSE 


noq+[401 


-1.19 ± [0.1] 


-1.83 ± [0.4] 


2 


25.00 ± 5.00 


40 


970828 


Konus 


900+ [30] 
^^°-[30] 


-0.704 ± [0.1] 


-2.07 ± [0.4] 


2 


2.20 ±0.40 


40 


971214 


BATSE 


190+ [20] 


-0.78 ± [0.1] 


-2.57 ± [0.4] 


2 






980703 


BATSE 


254+[25] 


-1.31 ± [0.1] 


-2.4 ± [0.4] 


2 


3.40 ± 0.50 


40 


990123 


BATSE 


eo4+[6o] 

°"^-[60] 


-0.9 ± [0.1] 


-2.48 ± [0.4] 


2 


2.04 ± 0.46 


41 


990506 


BATSE 




-1.37 ± [0.1] 


-2.15 ± [0.4] 


2 






990510 


BATSE 


126+''°' 


-1.28 ± [0.1] 


-2.67 ± [0.4] 


2 


1.60 ±0.20 


42 


990705 


Konus 


189115 


-1.05 ± 0.21 


-2.2 ± 0.1 


1 


1.00 ±0.20 


40 


991208 


Konus 


-^^"-[20] 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


3 






991216 


BATSE 


O1Q+130] 
"^^**-[30] 


-1.23 ± [0.1] 


-2.18 ± [0.4] 


2 


1.20 ±0.40 


40 


000131 


Konus 


163111 


-1.2 ±0.1 


-2.4 ±0.1 


4 






000210 


Konus 


408111 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


5 






000911 


Konus 


^^"-[lOO] 


-0.84 ± [0.1] 


[-2.2] ± [0.4] 


6 






000926 


Konus 


10017 


[-1.1] ± [0.4] 


-2.43 ± [0.4] 


5 






010222 


Konus 


309112 


-1.35 ±0.19 


-1.64 ±0.02 


5 


0.93 ±0.10 


40 


010921 


HETE 


89l?2 


-1.6 ±0.1 


[-2.2] ± [0.4] 


7 






020124 


HETE 


87ll« 


-0.81°:? 


[-2.2] ± [0.4] 


7 


3.00 ± 0.40 


43 


020405 


Konus 


3641^° 


[-1.1] ± [0.4] 


-1.87 ±0.2 


8 


1.67 ± 0.52 


44 


020813 


HETE 


14011^3 


-0.94 ± 0.03 


, r:7+0.03 


7 


0.43 ±0.06 


40 


021004 


HETE 




-1 ± 0.2 


[-2.2] ± [0.4] 


7 


4.74 ± 0.50 


45 


021211 


HETE 


461^ 


-0.9 ±0.1 


9 9+0.1 
^•^-0.3 


7 






030115 


HETE 


QQ+53 


-1.3 ±0.1 


[-2.2] ± [0.4] 


7 






030226 


HETE 


971?^ 


-0.9 ±0.2 


[-2.2] ± [0.4] 


7 


1.04 ±0.12 


46 


030323 


HETE 


441^^ 


-0.8 ± 0.8 


[-2.2] ± [0.4] 


9 






030328 


HETE 


13011^ 


-1.14 ±0.03 


9 1+0.2 
^■^-0.4 


7 


0.80 ±0.10 


47 


030329 


HETE 


681^ 


-1.26 ±0.02 


-2.28 ±0.06 


7 


0.50 ±0.10 


43 


030429 


HETE 


351^2 


_1 1+0-3 
^•^-0.2 


[-2.2] ± [0.4] 


7 


1.77 ±1.00 


48 


030528 


HETE 


32lg 


1 q+0.2 


9 7+0.3 
-1.0 


7 






040924 


HETE 


67l« 


-1.17 ± [0.1] 


[-2.2] ± [0.4] 


10, 11 






041006 


HETE 


63111 


-1.37 ± [0.1] 


[-2.2] ± [0.4] 


11 


0.16 ±0.04 


49 


050126 


Swift 


47lf 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


12 
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Table 2 — Continued 



GRB Satellite Ep^^k (keV)* a" Ref> tjet{day) Ref> 



050223 


Swift 


+ [101 
°^-[10] 


-1.46 ± [0.1] 


[-2.2] ± [0.4] 


13 






050315 


Swift 


39I7 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


14 






050401 


Swift 




-0.9 ±0.3 


-2.55 ± 0.3 


15 






050406 


Swift 


25l?| 


[-1.1] ± [0.4] 


-2.56 ± 0.35 


16 






050408 


HETE 


[100]+|^°°1 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 








050416A 


Swift 


^ --1-9 

15I3 


[-1.1] ± [0.4] 


-3.4 ±0.4 


17 






050505 


Swift 




-0.31 ± 1 


[-2.2] ± [0.4] 


12 






050525A 


Swift 


8i_i 


-1.01 ± 0.06 


-3.26 ± 0.2 


15 


0.40 ± 0.10 


50, 


050603 


Swift 


344lf2 


-1.03 ±0.06 


-2.03 ±0.1 


15 






050730 


Swift 


, -L9R 

124^26 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


14 






050802 


Swift 


-L9R 

121128 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


14 






050814 


Swift 


60l^* 


0±0.6 


[-2.2] ± [0.4] 


18 






050820A 


Swift 




-1.25 ±0.1 


[-2.2] ± [0.4] 


15 


188 ±2.00 


52 


050824 


Swift 


15t| 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


14 






050826 


Swift 


io5i:l^ 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


14 






050908 


Swift 


4ll| 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


12 






050922C 


Swift 


198122 


-0.95 ± 0.07 


[-2.2] ± [0.4] 


15 


0.11 ±0.03 


53 


051016B 


Swift 


^ .4-7 

241 7 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


14 






051022 


Konus 


„-l-99 

5101^2 


-1.18 ±0.02 


[-2.2] ± [0.4] 


19 


2.90 ±0.20 


54 


051109A 


Swift 


1 CI +130 
161-35 


-1 r)c + 0.27 
^•^°-0.36 


[—2.2] ± [0.4] 


20 


0.60 ± 0.10 


55 


051111 


Swift 


220™ 


-1±0.18 


[-2.2] ± [0.4] 


18 


0.0080 ± 0.0003 


56 


060108 


Swift 


65l«oo 


[-1.1] ±0.4 


[-2.2] ± [0.4] 


13 






060115 


Swift 


62118 


-1±0.3 


[-2.2] ± [0.4] 


21 






060206 


Swift 


78lf 


-1.2 ±0.18 


[-2.2] ± [0.4] 


22 


0.57 ± 0.06 


57, 


060210 


Swift 


149ir 


-1.18 ±0.03 


[-2.2] ± [0.4] 


13 






060223A 


Swift 


71+100 
'-'--10 


-1.18 ±0.31 


[-2.2] ± [0.4] 


13 






060418 


Swift 


230l[-| 


-1.5 ± [0.1] 


[-2.2] ± [0.4] 


23 






060502A 


Swift 


isel^" 


-1.18 ±0.15 


[-2.2] ± [0.4] 


13 






0605 lOB 


Swift 


Qr:+[60] 
^^-[30] 


-1.47 ±0.18 


[-2.2] ± [0.4] 


13 






060512 


Swift 


22l« 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


14 






060522 


Swift 


801^2^ 


-0.7 ±0.48 


[-2.2] ± [0.4] 


18 






060526 


Swift 




[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


13 


2.41 ± 0.06 


59 
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Table 2 — Continued 



GRB 


Satellite 


Bpeafc (keV)'' 


a'- 


pa. 




Re 


060604 


Swift 


40+151 


-1.34 ± [0.3] 


[-2.2] ± 


0.4] 


13 


060605 


Swift 




-0.3 ±0.42 


[-2.2] ± 


0.4] 


18 


060607A 


Swift 


120llf 


-1.06 ±0.18 


[-2.2] ± 


;o.4] 


13 


060707 


Swift 


mtf 


-0.6 ±0.42 


[-2.2] ± 


;o.4] 


22 


060714 


Swift 


103l?i 


[-1.1] ± [0.4] 


[-2.2] ± 


;o.4] 


24 


060729 


Swift 


Qltl 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


25 


060814 


Swift 


257lJt 


-1.43 ±0.09 


[-2.2] ± 


0.4] 


26 


060904B 


Swift 




-1±0.42 


[-2.2] ± 


0.4] 


18 


060908 


Swift 




-1 ±0.18 


[-2.2] ± 


;o.4] 


22 


060926 


Swift 


20111 


[-1.1] ± [0.4] 


[-2.2] ± 


;o.4] 


14 


060927 


Swift 


72tf 


-0.9 ±0.24 


[-2.2] ± 


;o.4] 


22 


061007 


Swift 


3991 


-0.7 ±0.02 


-2.61 ± 


0.09 


14 


061110A 


Swift 




[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


061110B 


Swift 


517153 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


061121 


Swift 


606155 


-1 or, + 0.02 

-"-•■^^-0.03 


[-2.2] ± 


0.4] 


27 


061222B 


Swift 


4911 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


070110 


Swift 


llOlfo 


[-1.1] ± [0.4] 


[-2.2] ± 


;o.4] 


28 


070208 


Swift 


5111S 


[-1.1] ± [0.4] 


[-2.2] ± 


;o.4] 


14 


070318 


Swift 


154l« 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


070411 


Swift 


83111 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


070506 


Swift 


311^ 


-5.00 ±1.82 


2.01 ±0.22 


18 


070508 


Swift 


233I7 


-0.96 ±0.13 


[-2.2] ± 


;o.4] 


14 


070521 


Swift 


222III 


-0.93 ±0.07 


[-2.2] ± 


;o.4] 


14 


070529 


Swift 


180l|2 


[-1.1] ± [0.4] 


[-2.2] ± 


;o.4] 


14 


070611 


Swift 


92112 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


070612A 


Swift 


8711? 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


070714B 


Swift 


11201^^3 


-0.86 ±0.06 


[-2.2] ± 


0.4] 


14 


070802 


Swift 


701^5 


[-1.1] ± [0.4] 


[-2.2] ± 


;o.4] 


14 


070810A 


Swift 


4412 


[-1.1] ± [0.4] 


[-2.2] ± 


;o.4] 


14 


071003 


Swift 


799l^f 


-0.97 ± 0.04 


[-2.2] ± 


;o.4] 


14 


071010A 


Swift 


97+10 


[-1.1] ± [0.4] 


[-2.2] ± 


0.4] 


14 


071010B 


Swift 


52l« 


, 9,.+0.4B 
^■^^-0.30 


-2 65+°'l® 


29 



tjet{day) Ref.' 



0.24 ± 0.02 60 







12 ± [0.01] 24 
50 ± 0.06 61 



1.16 ±[0.16] 62 



3.44 ± 0.39 63 
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Table 2 — Continued 



GRB Satellite Ep^ak (keV)'» (3^ Ref> tj^t{day) Ref.'' 

14 
30 
14 
14 
14 
31 
32 
14 
33 
34 
35 
14 
14 
14 
36 
37 
38 
14 
39 

''The values reported in square brackets are conservative estimations for uncertainties not reported 
in the original paper. 

''References. — ( 1 ) Amati et al. 2002; ( 2 ) Jimenez, Band, & Piran 2001; ( 3 ) Golenetskii 2005, 
private communication; ( 4 ) Andersen et al. 2000; ( 5 ) Ulanov et al. 2005; { 6 ) Price et al. 2002a; 
( 7 ) Sakamoto et al. 2005; ( 8 ) Price et al. 2003a; ( 9 ) Atteia et al. 2005; ( 10 ) Golenetskii et 
al. 2004; ( 11 ) HETE Bursts 2006, available at |http: / /space. mit .edu /HETE /Bursts / 1 ( 12 ) Krimm 
2005, private communication; ( 13 ) Krimm 2006, private communication; ( 14 ) calculated from the 
relation in Zhang et al. 2007; { 15 ) Krimm et al. 2006a; ( 16 ) Schady et al. 2006; { 17 ) Sakamoto 
et al. 2006; ( 18 ) Butler et al. 2007; ( 19 ) Golenetskii et al. 2005a; { 20 ) Golenetskii et al. 2005b; { 
21 ) Barbier et al. 2006a; ( 22 ) Sakamoto et al. 2008; ( 23 ) Golenetskii et al. 2006a; { 24 ) Krimm et 
al. 2007; ( 25 ) Rykoff et al. 2006; { 26 ) Golenetskii et al. 2006b; { 27 ) Golenetskii et al. 2006c; { 28 
) Amati et al. 2007; ( 29 ) Golenetskii et al. 2007a; ( 30 ) Golenetskii et al. 2007b; ( 31 ) Golenetskii 



071031 


Swift 


241^ 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


071117 


Swift 


_ IT Art 


-1.53 ± 0.09 


[-2.2] ± [0.4] 


071122 


Swift 


'"^-30 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080210 


Swift 


iyo + 15 

'"^-15 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080310 


Swift 




[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080319B 


Swift 


65lt| 


-0.82 ± 0.01 


~3-»'_0.66 


080319C 


Swift 


307l|^ 


-1.01 ± 0.08 


1 ov+0-09 


080330 


Swift 


20l| 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080411 


Swift 


259l?i 


1 CI +0.02 
^•3^-0.03 


[-2.2] ± [0.4] 


080413A 


Swift 


170l^« 


-1.20 ±0.06 


[-2.2] ± [0.4] 


080413B 


Swift 


'■^-10 


-1.26 ± 0.16 


[-2.2] ± [0.4] 


080430 


Swift 




[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080516 


Swift 




[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080520 


Swift 




[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080603B 


Swift 




-0 94+° '^3 
"•^^-0.45 


-1.96 ± [0.4] 


080605 


Swift 


246tlt 


-1.02 ± 0.06 


[-2.2] ± [0.4] 


080607 


Swift 


39411 


-1.06 ± 0.05 


[-2.2] ± [0.4] 


080707 


Swift 


yO+20 

''-'-20 


[-1.1] ± [0.4] 


[-2.2] ± [0.4] 


080721 


Swift 


485l« 


_0 93+0.06 
'^•^•^-0.05 


_2 43+0.15 
^•^'-'-0.25 
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et al. 2008a; ( 32 ) Golenetskii et al. 2008b; ( 33 ) Golenetskii et al. 2008c; ( 34 ) Ohno et al. 2008; 

( 35 ) Barthelmy et al. 2008; ( 36 ) Golenetskii et al. 2008d; ( 37 ) Golenetskii et al. 2008e; ( 38 ) 
Golenetskii et al. 2008f; ( 39 ) Golenetskii et al. 2008g; ( 40 ) Bloom et al. 2003b; { 41 ) Kulkarni et 
al. 1999; ( 42 ) Israel et al. 1999; ( 43 ) Berger et al. 2003; ( 44 ) Price et al. 2003a; ( 45 ) Holland 
et al. 2003; ( 46 ) Klose et al. 2004; ( 47 ) Andersen et al. 2003; ( 48 ) Jakobsson et al. 2004; ( 49 ) 
Stanek et al. 2005; ( 50 ) Liang et al. 2008; ( 51 ) Mirebel et al. 2007; ( 52 ) Cenko et al. 2006a; ( 53 
) Li et al. 2005; ( 54 ) Racusin et al. 2005; ( 55 ) Yost et al. 2007; ( 56 ) Guidozi et al. 2007; ( 57 ) 
Curran ot al. 2007; ( 58 ) Stcnik ct al. 2005; { 59 ) Thocnc ct al. 2008c; ( 60 ) Ferrero et al. 2009; ( 
61 ) Grupe et al. 2007; ( 62 ) Page et al. 2007; ( 63 ) Kann et al. 2007; 
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Table 3. Luminosity Relations 



Indicator 






Relation 




^6" 


/T * 

Vsys 


Tlag 


logL = 


52.21 - 


0.981og[n„g(l + z)-V0.1s] 


0.04 


0.03 


0.48 


V 


logL = 


51.06 + 


1.35 loglVCl + z)/0.02] 


0.09 


0.09 


0.87 


^peak 


logL = 


52.10 + 


1.881og[Epeafe(l + z)/300keV] 


0.03 


0.01 


0.57 




log E~f -- 


= 50.57 


+ 1.63log[Ep^ak{l + z)/300keV] 


0.09 


0.03 


0.16 




logL = 


52.84 - 


1.701og[ri{T(l + 2)-V0.1s] 


0.04 


0.05 


0.53 




log L > 


5n.:i2 + 


21<>s:..V,„„s] lor A-,„„t > 2 









cTa, CTjj are the 1 — a uncertainties of the fitting parameters, while a is the 
intercept on the y-axis, and 6 is the slope of the fitting line. (Tsys is the systematic 
uncertainty added in quadrature to the measurement errors, with which the X^ed °^ 
points about the best fit line is unity. 



Table 4. Comparison of Luminosity Relations 



Luminosity Relation 






pre-Swift burst 




Swift burst 


F 


^norm,dif f 




logL 


= 52.14 


-0.951og[Ti,9(l + z)-V0.1s] 


logL 


= 52.26 - 1.001og[n„g(l + 2)-V0.1s] 


0.98 


0.11 


Variability-L 


logL 


= 51.38 


+ 1.241og[V'(l + 2)/0.02] 


logL 


= 50.92 + 1.34 logfVCl + 2)/0.02] 


1.04 


0.09 




logL 


= 51.42 


+ 1.71 log[£'p,„fc(l + z)/300keV] 


logL 


= 51.13 + 1.88 log[L;pe^fe(l + z)/300keV] 


1.02 


0.08 




logL 


= 52.87 


-1.441og[THT(l + ^)-V0.1s] 


logL 


= 52.81 - 1.851og[TflT(l + ^)~VO-ls] 


1.00 


0.18 



^The <Tnorm,d%ff the 1 — (T range of the difference of normahzation values between Swift and pre-Swift luminosity relations. It is calculated 

yj ^norm,preSwift norm, Swift' 



I 

CO 

I 
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Table 5. Our Redshifts and Spectroscopic Redshifts. 



GRB 


Experiment 




Ref.'' 


Zb t 




Zmax 


N 


970228 


Konus 


0.84 


1 


0.38 


0.18 


0.86 


2 


970508 


BATSE 


0.84 


2 


1.51 


1.11 


2.09 


4 


970828 


Konus 


0.96 


3 


1.10 


0.73 


1.83 


4 


971214 


BATSE 


3.42 


4 


5.21 


2.00 


20 


3 


980703 


BATSE 


0.97 


5 


0.84 


0.66 


1.11 


4 


990123 


BATSE 


1.61 


6 


0.91 


0.67 


1.27 


4 


990506 


BATSE 


1.31 


7 


1.16 


0.77 


1.84 


4 


990510 


BATSE 


1.62 


8 


2.47 


1.69 


3.66 


5 


990705 


Konus 


0.84 


9 


0.63 


0.45 


0.94 


4 


991208 


Konus 


0.71 


10 


0.46 


0.29 


0.78 


3 


991216 


BATSE 


1.02 


11 


0.60 


0.47 


0.79 


5 


000131 


Konus 


4.50 


12 


1.75 


0.90 


4.41 


3 


000210 


Konus 


0.85 


13 


0.59 


0.35 


1.04 


2 


000911 


Konus 


1.06 


14 


2.88 


1.06 


18.31 


3 


000926 


Konus 


2.07 


15 


0.60 


0.31 


1.41 


2 


010222 


Konus 


1.48 


16 


0.45 


0.35 


0.61 


4 


010921 


HETE 


0.45 


17 


0.40 


0.28 


0.59 


3 


020124 


HETE 


3.20 


18 


5.27 


3.24 


8.12 


5 


020405 


Konus 


0.70 


19 


0.90 


0.64 


1.32 


4 


020813 


HETE 


1.25 


20 


0.98 


0.72 


1.36 


5 


021004 


HETE 


2.32 


21 


2.41 


1.53 


4.06 


5 


021211 


HETE 


1.01 


22 


0.71 


0.49 


1.07 


3 


030115 


HETE 


2.50 


23 


2.24 


1.34 


4.12 


3 


030226 


HETE 


1.98 


24 


2.39 


1.59 


3.58 


5 


030323 


HETE 


3.37 


25 


1.52 


0.55 


20 


2 


030328 


HETE 


1.52 


26 


1.28 


0.80 


2.34 


5 


030329 


HETE 


0.17 


27 


0.24 


0.19 


0.29 


5 


030429 


HETE 


2.66 


28 


4.62 


2.41 


8.25 


4 


030528 


HETE 


0.78 


29 


0.77 


0.53 


1.14 


3 


040924 


HETE 


0.86 


30 


0.61 


0.42 


0.91 


3 


041006 


HETE 


0.71 


31 


0.64 


0.48 


0.89 


4 


050126 


Swift 


1.29 


32 


1.64 


1.00 


2.91 


3 
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Table 5 — Continued 



GRB Experiment Zspec Ref.* Z(,est Zmin Zmax Nrel 



050223 


Swift 


0.59 


33 


2.80 


0.94 


20 


1 


050315 


Swift 


1.95 


34 


0.82 


0.45 


1.78 


3 


050401 


Swift 


2.90 


35 


2.02 


1.23 


3.63 


4 


050406 


Swift 


2.44 


36 


1.68 


0.71 


20 


2 


050408 


HETE 


1.24 


37 


1.43 


0.90 


2.40 


3 


050416A 


Swift 


0.65 


38 


0.73 


0.43 


1.33 


2 


050505 


Swift 


4.27 


39 


2.73 


1.68 


4.67 


4 


050525A 


Swift 


0.61 


40 


0.92 


0.71 


1.21 


5 


050603 


Swift 


2.82 


41 


1.24 


0.69 


2.69 


2 


050730 


Swift 


3.97 


42 


> 2.53 






2 


050802 


Swift 


1.71 


43 


1.09 


0.61 


2.41 


3 


050814 


Swift 


5.30 


44 


2.70 


0.89 


20.00 


2 


050820A 


Swift 


2.61 


45 


2.97 


1.96 


4.80 


4 


050824 


Swift 


0.83 


46 


0.61 


0.30 


1.50 


1 


050826 


Swift 


0.30 


47 


> 2.42 






2 


050908 


Swift 


3.35 


48 


2.02 


0.84 


9.23 


2 


050922C 


Swift 


2.20 


49 


1.51 


1.20 


1.90 


5 


051016B 


Swift 


0.94 


50 


0.65 


0.32 


1.52 


2 


051022 


Konus 


0.80 


51 


0.95 


0.66 


1.51 


3 


051109A 


Swift 


2.35 


52 


1.33 


0.74 


2.34 


3 


051111 


Swift 


1.55 


53 


1.09 


0.79 


1.48 


4 


060108 


Swift 


2.03 


54 


2.96 


0.74 


20 


2 


060115 


Swift 


3.53 


55 


2.63 


0.99 


20 


3 


060206 


Swift 


4.05 


56 


1.73 


0.93 


4.02 


4 


060210 


Swift 


3.91 


57 


2.54 


1.95 


3.97 


4 


060223A 


Swift 


4.41 


58 


> 2.06 






3 


060418 


Swift 


1.49 


59 


1.41 


0.92 


2.25 


4 


060502A 


Swift 


1.51 


60 


0.86 


0.55 


1.37 


3 


060510B 


Swift 


4.90 


61 


> 1.72 






2 


060512 


Swift 


0.44 


62 


0.64 


0.31 


1.69 


1 


060522 


Swift 


5.11 


63 


7.24 


1.08 


20 


1 


060526 


Swift 


3.21 


64 


3.50 


2.07 


6.57 


5 
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Table 5 — Continued 



GRB Experiment Zspec Ref.'' Zbest ^min Zmax Nrel 



060604 


Swift 


2.68 


65 


2.72 


0.97 


20 


2 


060605 


Swift 


3.80 


66 


6.79 


1.96 


20 


4 


060607A 


Swift 


3.08 


67 


1.47 


0.92 


2.48 


3 


060707 


Swift 


3.43 


68 


2.15 


0.78 


20 


2 


060714 


Swift 


2.71 


69 


3.77 


1.09 


20 


3 


060729 


Swift 


0.54 


70 


1.66 


0.81 


3.48 


3 


060814 


Swift 


0.84 


71 


0.93 


0.63 


1.43 


4 


060904B 


Swift 


0.70 


72 


1.89 


0.81 


5.46 


3 


060908 


Swift 


2.43 


73 


2.13 


1.30 


3.84 


4 


060926 


Swift 


3.21 


74 


1.44 


0.83 


2.71 


3 


060927 


Swift 


5.60 


75 


3.24 


1.81 


6.44 


4 


061007 


Swift 


1.26 


76 


0.94 


0.65 


1.42 


4 


061110A 


Swift 


0.76 


77 


> 1.82 






1 


061110B 


Swift 


3.44 


78 


8.28 


3.35 


20 


4 


061121 


Swift 


1.31 


79 


0.79 


0.58 


1.11 


5 


061222B 


Swift 


3.36 


80 


1.10 


0.50 


3.68 


2 


070110 


Swift 


2.35 


81 


> 1.995 






1 


070208 


Swift 


1.17 


82 


1.76 


0.70 


10.15 


2 


070318 


Swift 


0.84 


83 


3.42 


1.43 


18.32 


2 


070411 


Swift 


2.95 


84 


2.73 


1.53 


4.72 


2 


070506 


Swift 


2.31 


85 


2.73 


1.80 


4.38 


3 


070508 


Swift 


0.82 


86 


1.25 


0.81 


2.02 


4 


070521 


Swift 


0.55 


87 


2.08 


1.27 


3.73 


4 


070529 


Swift 


2.50 


88 


> 1.67 






1 


070611 


Swift 


2.04 


89 


5.55 


1.24 


20 


1 


070612A 


Swift 


0.62 


90 


1.82 


0.85 


5.61 


3 


070714B 


Swift 


0.92 


91 


3.12 


1.74 


6.97 


3 


070802 


Swift 


2.45 


92 


9.37 


1.51 


20 


1 


070810A 


Swift 


2.17 


93 


1.53 


0.97 


2.56 


3 


071003 


Swift 


1.10 


94 


0.83 


0.56 


1.25 


4 


071010A 


Swift 


0.98 


95 


0.86 


0.38 


2.88 


1 


071010B 


Swift 


0.95 


96 


0.76 


0.53 


1.12 


4 
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Table 5 — Continued 



GRB 


Experiment 


Zspec 


Ref.^ 






Zmax 


Nrel 


071031 


Swift 


2.69 


97 


1.10 


0.52 


3.40 


2 


071117 


Swift 


1.33 


98 


1.13 


0.76 


1.74 


3 


071122 


Swift 


1.14 


99 


> 1.05 






1 


080210 


Swift 


2.64 


100 


2.35 


1.37 


4.42 


4 


080310 


Swift 


2.43 


101 


1.19 


0.60 


3.22 


3 


080319B 


Swift 


0.94 


102 


1.45 


0.95 


2.34 


4 


080319C 


Swift 


1.95 


103 


3.62 


1.31 


20 


3 


080330 


Swift 


1.51 


104 


0.79 


0.37 


1.94 


2 


080411 


Swift 


1.03 


105 


0.51 


0.36 


0.74 


4 


080413A 


Swift 


2.43 


106 


2.42 


1.47 


4.36 


4 


080413B 


Swift 


1.10 


107 


0.73 


0.50 


1.11 


3 


080430 


Swift 


0.77 


108 


1.56 


1.00 


2.60 


3 


080516 


Swift 


3.20 


109 


3.62 


1.86 


7.96 


3 


080520 


Swift 


1.55 


110 


0.52 


0.25 


1.21 


1 


080603B 


Swift 


2.69 


111 


4.55 


2.40 


9.88 


4 


080605 


Swift 


1.64 


112 


1.06 


0.71 


1.65 


4 


080607 


Swift 


3.04 


113 


1.32 


0.82 


2.29 


4 


080707 


Swift 


1.23 


114 


2.35 


0.77 


20 


2 


080721 


Swift 


2.60 


115 


1.13 


0.72 


1.86 


4 
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2006; ( 24 ) Price et al. 2003b; ( 25 ) Vreeswijk et al. 2003; ( 26 ) Martini et al. 
2003; ( 27 ) Greiner et al. 2003; ( 28 ) Weidinger et al. 2003; ( 29 ) Ran, Salvato, 
& Greiner 2005; ( 30 ) Wiersema et al. 2004; ( 31 ) Fugazza et al. 2004; ( 32 ) 
Berger et al. 2005a; ( 33 ) Berger & Min-Su Shin 2006; ( 34 ) Kelson & Berger 



-54- 



2005; ( 35 ) Fynbo et al. 2005a; ( 36 ) Schady et al. 2006; ( 37 ) Berger et al. 
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Piranomonte et al. 2006; ( 56 ) Fynbo et al. 2006a; ( 57 ) Cucchiara et al. 2006a; 
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( 61 ) Price et al. 2006; ( 62 ) Bloom et al. 2006a; ( 63 ) Ccnko ct al. 2006b; ( 64 
) Berger & Gladders 2006; ( 65 ) Castro-Tirado et al. 2006; ( 66 ) Peterson et al. 
2006; ( 67 ) Ledoux et al. 2006; ( 68 ) Jakobsson et al. 2006a; ( 69 ) Jakobsson et 
al. 2006b; ( 70 ) Thoene et al. 2006a; ( 72 ) Fugazza et al. 2006; ( 73 ) Rol et al. 
2006; ( 74 ) D'Elia, et al. 2006; ( 75 ) Fynbo et al. 2006b; ( 76 ) Osip et al. 2006; 
( 77 ) Thoene et al. 2006b; ( 78 ) Fynbo ct al. 2006c; ( 79 ) Bloom et al. 2006; ( 
80 ) Berger 2006; ( 81 ) Jaunsen et al. 2007a; ( 82 ) Cucchiara et al. 2007; ( 83 ) 
Jaunsen et al. 2007b; ( 84 ) Jakobsson et al. 2007a; ( 85 ) Thoene et al. 2007a; ( 
86 ) Jakobsson et al. 2007b; ( 87 ) Hattori et al. 2007; ( 88 ) Berger et al. 2007; ( 
89 ) Thoene et al. 2007b; ( 90 ) Cenko et al. 2007a; ( 91 ) Graham et al. 2007; ( 
92 ) Prochaska et al. 2007a; ( 93 ) Thoene et al. 2007c; ( 94 ) Perley et al. 2007; 
( 95 ) Prochaska, J. X. et al. 2007b; ( 96 ) Ccnko ct al. 2007b; ( 97 ) Ledoux 
et al. 2007; ( 98 ) Jakobsson ct al. 2007c: ( 99 ) Cucchiara et al. 2007b; ( 100 ) 
Jakobsson et al. 2008a; ( 101 ) Prochaska, et al. 2008a; ( 102 ) Vreeswijk et al. 
2008a; ( 103 ) Wiersema et al. 2008; ( 104 ) Malesani et al. 2008; ( 105 ) Thoene 
et al. 2008a; ( 106 ) Thoene et al. 2008b; ( 107 ) Vreeswijk et al. 2008b; ( 108 ) 
Cucchiara & Fox 2008; ( 109 ) Filgas et al. 2008; ( 110 ) Jakobsson et al. 2008b; 
( 111 ) Fynbo et al. 2008a; { 112 ) Jakobsson ct al. 2008c; ( 113 ) Prochaska, et 
al. 2008b; ( 114 ) Fynbo et al. 2008b; ( 115 ) D'Avanzo et al. 2008; 
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Fig. 1. — Tiag calculation - Cross-Correlation versus Offset plots. The left panel is the plot 
for GRB051111, which is a moderately bright burst without much noise in the light curve. 
But we can see that although the curve is relatively smooth, the peak region is quite 'broad' 
that we still need an accurate fitting of the curve, especially for the central region of the 
curve. Had we simply taken the offset of the highest measured cross-correlation, we would 
have the peak at -14 bins (with Tiag equal to 0.90 seconds), with this value being determined 
by noise in the cross correlation. Instead, we fitted a 7-order polynomial (smooth curve) and 
found a peak at -23.1 bins (with Tiag equal to 1.48 seconds), and this model fit is avoiding 
most of the random noise in the cross correlation. Note, for our choice of bands, the lag will 
be negative in bins but expressed as positive in time. The right panel is for GRB050820, 
which is a noisy burst. There is much scatter in the cross-correlation versus offset plot, 
from which we cannot determine the position of the peak without making a fit around the 
peak region. In this case, the highest cross correlation is at +23 bins (corresponding to a 
negative lag of -1.47 seconds), with the model fit (a 7-order polynomial also) providing a 
more reasonable peak at -15.7 bins {riag — 1.00 seconds). 
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Fig. 2. — The calculation of t^t and Np^ak- We can see an apparent 'plateau' in the lower 
plot, for box smoothing widths of more than 5. By simply extrapolating the plateau back to 
the y-axis (see dashed line), we get a minimum rise time (see circle at intercept on y-axis) 
which is appropriate for zero smoothing, and which is not affected by the Poisson noise. 
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Fig. 3. — Luminosity relations for long GRBs. riag, V, trt and Np^ak values are generated 
from light curves of the GRBs, while -Bpeafc values are generated from the spectra. The 
1 — a measurement uncertainties are used as the error bar. The best fit relation of the six 
is the Epeak — relation (Ghirlanda's relation), while the most scattered one is the V — L 
relation. However, with the restrictions on the measurements of Epeak (especially for Swift 
bursts) and jet break time, not many bursts can be included in the Ghirlanda's relation. 
The Npeak — L relation provides only a lower limit on the luminosity, but this is a strong 
hmit never violated, for example where a 10-peak burst always has L > 10^^-^^ ergs. 
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Fig. 4. — Comparison of the luminosity relations for pre-Swift and Swift bursts. The filled 
diamonds are for pre-Swift bursts, while the open squares are for Swift bursts. From the plot 
we can tell that the luminosity relations are identical between pre-Swift and Swift bursts to 
within the usual uncertainties. Ghirlanda's relation is not compared because the sample of 
bursts with both Ep^ak and jet break time is too small. 
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Fig. 5. — One example of the plot in our calculation. Taking the ^cos as model value and 
each of the Hi as observed value, we can have the fo^^ each of the indicators, as shown in 
the left panel. Each of the luminosity relations corresponding with these curves are labeled 
just beside the lines, and the vertical line on the left is the luminosity lower limit given by 
the Npeak — L relation. We sum up all the included x^s to get a total plot, as shown in 
the right panel. Then the minimum corresponds with our Zhest, and = Xmin + 1) 4, 9 
gives us the l,2,3-(7 range of our estimated redshifts. 
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log(z^c) 

Fig. 6. — Comparison of our redshifts with spectroscopic redshifts. The diamonds are our 
Zbest and the error bars gives the 1-a redshift range. We can see that 70% of the spectroscopic 
redshifts fall in our 1-a range, and our Zbest scatter uniformly around the spectroscopic 
redshift (with {log {zbest / z spec)) = 0.01). This demonstrates that our method has accurate 
error bars and negligible biases. 
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log(Zspec) 

Fig. 7. — We made a comparison between our Zbest and Zbestji- Zbest,ii are the calculated 
redshift values when we apply Swift luminosity relations on pre- Swift data and pice- Swift 
relations on Swift data. The analysis on Zspecji and Zspec data shows that the average value 
of logio{zbest,ii/ Zspec) IS "0.02, and the RMS scatter of logio{zbest,ii / Zspec) is 0.27, both of 
which are equal to those of our z^est value within error bars. From this we see that these 
two sets of redshifts do not differ much from each other, which also means that the pre-Swift 
luminosity relations and Swift luminosity relations do not have significant difference. It also 
tells us that the effect of redshift involved in our calculation is negligible. 



